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We investigate effects of cluster particle correlations on two- and three-particle azimuth correlator 
observables sensitive to local strong parity violation. We use two-particle angular correlation mea- 
surements as input and estimate the magnitudes of the effects with straightforward assumptions. 
We found that the measurements of the azimuth correlator observables by the STAR experiment 
can be entirely accounted for by cluster particle correlations together with a reasonable range of 
cluster anisotropy in non-peripheral collisions. Our result suggests that new physics, such as local 
strong parity violation, may not be required to explain the correlator data. 
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I. INTRODUCTION 

Relativistic heavy-ion collisions at RHIC have created a hot and dense medium that exhibits properties of a strongly 
coupled Quark-Gluon Plasma (sQGP) Approximate chiral symmetry may be restored in the sQGP. It was recently 
suggested that metastable domains may form in such sQGP state where the parity and time-reversal symmetries are 
locally violated 0, [j| . Such violation would lead to separation of positive and negative particles due to the chiral 
magnetic effect along the system's orbital angular momentum into the two hemispheres separated by the reaction 
plane [1, i. 

The most direct consequence of this charge separation is a negative correlation of multiplicity asymmetry of positive 
particles between the hemispheres separated by the reaction plane and that of negative particles [2J . Such a negative 
correlation is in addition to any (background) correlations that may exist due to other dynamics of the collision, the 
magnitude of which may be assessed by correlation of multiplicity asymmetries between the hemispheres separated 
by the plane normal to the reaction plane. 

Another consequence of the charge separation is a positive correlator (cos(ev + (3 — 2^)) of unlike-sign (US) particle 
pairs and a negative (cos(a + f3 — 2ip)} of like-sign (LS) particle pairs, where a and /3 are the azimuthal angles of 
the two particles and i\) is the reaction plane azimuth |4[. The reaction plane azimuthal angle is, however, not fixed 
but random in heavy-ion collisions. In order to estimate the reaction plane angle, a third particle, c, may be used to 
correlate with a and /3, correcting for the resolution effect (i>2, C j the elliptic flow of particle c) Namely 

(cos(a + 0- 2if>)) « (cos(a + j3 - 2c))/v 2 ,c- (1) 

This assumes three-particle correlation is negligible. 

The three-particle azimuthal correlator (cos(a + (3 — 2c)) has been measured and is used to deduce the two-particle 
azimuthal correlator (cos(a + (3 — 2ij))) by the STAR experiment 0. The measurements show a negative correlator 
(cos(a + (3 — 2ip)) for LS pairs and a small, close to zero, correlator for US pairs [5j. The LS pair result is qualitatively 
consistent with the expectation from local strong parity violation @, The US pair result, however, is inconsistent 
with the initial expectation where the US and LS pair correlations should be equal in magnitude and opposite in 
sign [2| . To explain the preliminary version of the STAR data [6| , it was suggested that at least one of the particles 
from a back-to-back US pair from local parity violation would have to traverse and interact with the medium and 
its angular correlation with the other particle of the pair would be significantly reduced In fact, in this medium 
interaction scenario, correlations between particle pairs from local parity violation domains formed in the interior of 
the collision medium would be lost, and only those from pairs emitted from the surface could survive. In other words, 
the LS pair correlation is due to those pairs from local parity violation domains on the surface, and the back-to-back 
US pair correlation is lost 

The three-particle correlator observable is parity-even and is subject to background correlations that are reaction- 
plane-dependent, some of which are discussed in detail in Ref. Q- This is easy to see in the following extreme: a 
small opening angle pair perpendicular to the reaction plane is indistinguishable from a back-to-back pair parallel to 
the reaction plane in the correlator variable cos(a + f3 — 2-0), and vice versa. More modestly, particle correlations 
from clusters which themselves possess anisotropy can give rise to observable signals in (cos(cv + (3 — 2ip)). Cluster 
correlations can have different effects on LS and US pairs, because LS and US contributions from clusters are likely 
different, such as from jet-correlations 0j- 

In this paper, we investigate quantitatively effects of cluster particle correlations on the azimuth correlator ob- 
servables (cos(a + (3 — 2c)) and (cos(a + (3 — 2i/j)). We first present analytical results. We then use experimental 
measurements of two-particle angular correlations H| as input to estimate effects of cluster particle correlations 
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on the correlator observables. Since we do not have experimental information on cluster anisotropy, we calculate 
how much cluster anisotropies are needed in order to fully account for the correlator measurements We then 
judge the plausibility of the needed cluster anisotropies to either confirm or disprove cluster correlations as a possible 
explanation for the correlator measurements. 



II. THE CLUSTER MODEL AND RESULTS 



Assume events are composed of hydrodynamic-like particles plus small-angle (SA, \Acf>\ < tt/2) clusters and back- 
to-back (BB, |A0| > 7r/2) clusters. We will use operational definitions of SA-cluster to be composed of one or more 
small-angle particle pairs, and BB-cluster to be composed of one or more back-to-back particle pairs. Note, with 
these operational definitions, that a conventional back-to-back cluster of a particles on one side and b particles on the 
other side is made of two SA-clusters (with numbers of particle pairs of a 2 and b 2 ) and one BB-cluster (with number 
of particles 2ab). Here we have taken a and b to be large for simplicity, and we will assume Poisson statistics for 
particle multiplicity in clusters. Also note that not all SA-clusters have a back-side partner. The anisotropies of SA- 
and BB-clusters can therefore be different. 

Now consider US and LS particle pairs from clusters. They can come from either SA-clusters (i.e. SA particle pairs) 
or BB-clusters (i.e. BB particle pairs) . The relative fractions of particle pairs from SA-clusters and BB-clusters can 
be different for US-pairs and LS-pairs. Let x us be the fraction of US-pairs from SA-clusters (and 1 — x us the fraction 
from BB-clusters), and x LS be the fraction of US-pairs from SA-clusters (and 1 — x LS the fraction from BB-clusters). 
We shall first obtain x us and x LS from two-particle angular correlation measurements. 

The following two-particle correlators are measured for US and LS pairs 

(cos(a- P)) vs = x us (cos(a- (3)) SA + (1 - x us )(cos(a - f3)) BB = x us w SA - (I - x us )w BB , (2) 

(cos(a - p)) LB = x LS (cos(a~ (3)) SA + (I - x LS )(cos(a - f3)) BB = x LS w SA - (f - x LS )w BB , (3) 

where 

w SA = (cos(a-/9)) SA , (4) 

w BB = -(cos(a-/3)) SB (5) 

are the average angular spread of particle pairs from SA- and BB-clusters, respectively. In Eqs. ((2]) and ([3]) we have 
taken the two-particle back-to-back correlations to be the same between US and LS pairs. We have assumed in 
Eqs. ([2]) and ((3j) that the SA two-particle azimuthal correlations to be of the same shape for US and LS pairs. This 
is a reasonable assumption because same-side correlations of US and LS pairs have similar shapes although their 
magnitudes are different, e.g. in jet-like correlations Q- 

Information about cluster particle pairs can be obtained from two-particle azimuthal correlations. STAR has 
measured two-particle correlations integrated over transverse momentum (pt), in (^A)</>a), the two-particle pseudo- 
rapidity and azimuth differences Q . The correlation functions are parameterized by the sum of a near-side Gaussian, 
a negative dipole, and a quadrupole corresponding to elliptic flow [8(. The sum of the two former terms is considered 
to be correlations due to clusters. It is given by 8]: 

d2N V ° ( <&\n( ^ A a (r\ 

ex P _ ^TT G (VA) - cos^a- (6) 



#A*?A \/2~K<7 V 2cr2 

Here G(?7a) is a Gaussian in 77A normalized to unity that is of no interest in our study. The first term in Eq. © 
r.h.s. is the near-side Gaussian and the second term is the negative dipole. We can obtain the SA pair azimuthal 
spread as 
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where A(t]a) is the two-particle tja acceptance of the STAR detector, and V is the integrated volume within the 
acceptance (which is not equal to Vb) 0,0- The extracted w SA ranges from 0.58 in peripheral Au+Au collisions to 
0.85 in central collisions. For w BB we take the width of the away-side dipole: 

§T/2 COS <M) COS <M#A 7T 

w BB =-{cos<Pa) bb = -3^- — = -■ (8) 

J tt/2 A M cos <p A d(p A * 
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The w BB is a fixed value because the away-side correlation shape can be satisfactorily described by the same functional 
form of a negative dipole. 

It is worthwhile to note that the cluster shape quantities are extracted from the measured angular correlations, thus 
they are immune to the underlying physics mechanisms generating those correlations. Besides cluster correlations, a 
negative dipole can be also generated, for example, by the statistical global momentum conservation. It was estimated, 
however, that the global momentum conservation effect is significantly smaller than the measured dipole strength. 
On the other hand, in order to estimate the cluster size that will be needed for our study below, some production 
mechanisms for the clusters have to be assumed. We will discuss those assumptions later. 

The pair quantities in Eqs. (J2)) and (J3j) measured in experiment are diluted by hydro-like particles (those pair 
quantities are all zero for hydro-like particle pairs as well as cross-pairs of hydro-like and cluster particles) : 



(cos(a - 
(cos(a - 



\ mcas 



f2,U S ( C0S ( a - 

L LS {cos(a- 



■P))vb, 
P)) LS - 



(9) 
(10) 



The dilution factors are 



N us /(N 2 /2), 
N LS /(N 2 /2) 7 



(11) 
(12) 



for US and LS pairs (numbers of pairs N us and N LS ), respectively, where we have assumed the total numbers of US 
and LS pairs are equal in the event (N is total particle multiplicity). The total number of cluster particle pairs is 



N US +N LS = N clvLSt M 2 lnst , 



(13) 



where iV c i ust is the number of clusters and -Mclust is the particle multiplicity per cluster (cluster size) . The numbers of 
US and LS pairs from clusters are not necessarily equal. Since there is no charge-sign difference in the back-to-back 
particle pair correlations: 



(1 x us)^v 



(1 



(14) 



we obtain the dilution factors: 



1 - ( X US +%L S )/ 2 

1 ~ x us 

1 - ( X US + X Ls)/ 2 



/a, 



/a, 



(15) 
(16) 



with 



h = N clust M 2 lust /N 2 



(17) 



Note N c i ust M 2 luEt /N -the number of correlated pairs per charged particle- is a measured quantity in two-particle 
correlation function; it is obtained by integrating the correlation function over the measured acceptance (which is the 
V in Eq. (J6]), V = N c \ ust M 2 lust /N) [1, H|. The measured V ranges from 0.3 in peripheral Au+Au collisions to 1.9 in 
central collisions. The measured N c i ust M 2 lust is not affected by the individual quantities of the cluster size M c i ust or 



the number of clusters iV c i UB t which are not experimentally measured 
We can now obtain x us and x LS from the measured (cos(a — /?)) 
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Eqs. © and © become 
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(18) 
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(20) 
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By simple algebra, we have 



wj A +w SA (3t us -t LS )/2 + 2w 



SA ^BB 



tV'us t LS ) 
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SA^US W Bb(^US t LS )/2 



(w SA + w BB )(w SA + (t us - t LS )/2) 
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0, 

(22) 
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x 2 _ w2 SA + w sa{^ls ~ tus)/ 2 + 2 W SA W BB + W BB (t LS -t u8 ) x W SA W BB +W SA t LS +W BB (t LS -t us )/2 = Q 

(w SA +w BB )(w SA + {t LS -t us )/2) Xls {w SA +w BB )(w SA +(t LS -t us )/2) 

(23) 

From Eqs. (|22t and (|23|) we can solve for x us and x LS . 

Table U shows the obtained fractions x us and x LS given cluster particle correlation inputs from two-particle angular 
correlations Q and two-particle corr elator quantities (cos(a - /3))™° as and (cos(a - /?))™ oas measured by STAR @. 

A. Effect of three-particle correlation from clusters 

STAR has measured three-particle correlators for US and LS pairs, a and /?, with a third particle c regardless of 
its charge sign. It is assumed that three-particle correlation is negligible, so particle c can be used as a single-particle 
estimator of the reaction plane to obtain the two-particle correlators from the three-particle correlator measurements 
by Eq. (TTJ). One supporting evidence for the assumption comes from the consistent results of (cos(a + j3 — 2tp)) rj 
(cos(a + f3 — 2c))/u2,e using particle c from the main Time Projection Chamber (TPC) or the forward TPC's while 
the particle pairs (US and LS) a and (3 come from the main TPC §. However, it is possible that the probability for 
a triplet to be correlated may drop with the pseudorapidity gap between the particle c and the other two particles in 
the main TPC, in a similar way to the ^2,c dependence on pseudorapidity [10j. 

If large clusters exist, as suggested by \ow-px two-particle angular correlation measurements Q, then finite three- 
particle correlation should exist. We shall estimate the effect of three-particle correlation using two-particle correla- 
tions Consider particle triplets from the same cluster, where SA and BB still stand for a pair of a and f3, and the 
third particle c can be on either side. 

(cos(a + p- 2c)) SA = (cos(A0 Q + Afa + 20 clust - 2c)) SA = (cos(A0 Q + Atfy? + 2A0 C )) SA , (24) 
(cos(a + p- 2c)) BB = (cos( A0 Q + Afo + 20 clust - 2c)) BB = (cos(A0 Q + Acf> + 2A0 C )) BB . (25) 



Here Ac/) = cf> — <^> c iust is the particle azimuth relative to the cluster axis, </> c iust- In Eq. (|24[) . which side the particle 
c is does not really matter (because the angle is 2c). In Eq. (j2"5j) . the particle c is on the same side of either a or 
P, and thereby either A<f> a or Aifip is larger than tt/2. Assuming that emission of daughter particles in clusters is 
independent of each other, we can obtain 

(cos(a + P - 2c)) SA « (cos Acj) a ) 2 sA (cos 2A4> C ) SA « (cos(A0 Q - Acf> )) SA (cos 2A4> C ) SA = w SA w 2 SA , (26) 
(cos(a + ^- 2c)) BB « (cosA0 Q ) BB (cos2A0 c ) gA « (cos(A0„ - A0 /3 )) BB (cos2A0 c ) SA = -w BB w 2SA . (27) 

Here A<f> a — A<fip = (j) A is the azimuth difference used in two-particle correlation measurement, and w 2 SA = 
(cos2A</>) SA is the average azimuthal spread of SA-clusters: 



Wn SA = (cos2A0) SA « (cos 20 



A 



,1/2 



1/2 

-^^ ex p{-^) G ( r lA)- A M coscj) A cos2(j} A dcj) A A(ri A )dr] A \ ' (yer^" 1 - 4A 0A /3\ 1/2 

(28) 



[^exp(-^)G(r, A )- Accost A ] dct> A A(ri A )dri A J \ v ~ 

Note in Eq. (|27p. it is the azimuthal spread of SA (not BB) clusters as well because particle c is always on the same 
side of either particle a or p. 

We can estimate effects of three-particle correlations in US and LS pairs of particles a and p by 



(cos(a + P-2c)) us = x us (cos(a + p-2c)) SA + (1 - x us )(cos(a + P - 2c)) BB x us w SA w 2SA - (l-x us )w BB i 
(cos(a +P- 2c)) LS = x LS (cos(a + P~ 2c)) SA + (1 - x LS )(cos(a + P~ 2c)) BB x LS w SA w 2SA - (1 - x LS )w BB w 2 0) 

Comparing Eqs. ([29]) and ([30j to Eqs. © and ©, we see that 

(cos(a + ^-2c)) [7S = (cos(a- P)) us -w 2 SA , (31) 
(cos(a + /3- 2c)) LS = (cos(a-P)) LS ■ w 2SA . (32) 

The only assumption in arriving at Eqs. pip and (|32"1) is that particle emission azimuths within a cluster are indepen- 
dent of each other. This is a reasonable assumption when the clusters consist of a relatively large number of particles. 
Under this assumption, three-particle correlation is completely determined by two-particle correlation. 
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FIG. 1: (Color online) The measured three-particle correlators for US and LS particle pairs (open data points) (fj, the estimated 
three-particle correlation effects (thin lines) using inputs from measurements of two-particle angular correlations ^ and two- 
particle correlators [H, and the remaining three-particle correlator magnitudes (solid data points), i.e. difference between open 
points and the thin lines. The left panel shows the three-particle correlators themselves versus centrality bin, and the right 
panel shows the number of participants (-/V par t) scaled three-particle correlators versus -/Vp ar t . 



These three-particle correlation effects are diluted by hydro-like particles, 

( C0S (a + /3-2c))^ as = f^ us (cos(a + (3~2c)) us , (33) 

(cos(a + /3-2c))™ as = / 3 ,Ls(cos(a + /3-2c)) iS , (34) 

by a factor of 

f3 ' US = l-(I~ + ^ s )/2 /3 ' (35) 

f3 ' LS = l"(^X)/2 /3 ' (36) 

with 

h = N clust M* lust /N 3 . (37) 

Again, the measured three-particle correlation is determined by the measured two-particle correlation: 

(cos(a + 0- 2c))™ as = (cos(a - /?))™ as • w 2 , SA ■ M clust /N, (38) 
(cos(a + /3- 2c))^» = (cos(a - /3))™ oas • w 2<SA ■ M clust /N. (39) 

The only inputs for this determination are the near-side angular spread w 2 SA by Eq. (f2"5f which is well measured and 
the cluster size M clust which can be estimated. We estimate M clust from the measured N c i ust Al^ lust assuming binary 
scaling for the number of clusters, N c \ ust . 

Table U shows the effect in three-particle azimuthal correlator estimated from the measured two-particle azimuthal 
correlator by Eqs. (|3"8")l and (|3T)]) . using cluster inputs from two-particle correlation measurements Figure [1] shows 
the estimated three-particle correlation effects in thin lines together with the measured three-particle correlator data 
in open points [5(. The estimated three-particle correlation effects from clusters are significantly larger than the US 
measurement (cos(a + f3 — 2c))™g as [5|. This implies that there must be some cancellation of US correlation in the 
data from other effects (one candidate is two-particle correlation as we will discuss below). Those estimated for LS 
pairs are smaller than measurement (cos(a + (3 — 2c))'°° as [5| in most of the centralities, by a factor of a few. 

The three-particle correlation effects should be first removed from the three-particle correlator measurements. After 
removing the estimated three-particle effects, the three-particle correlators (from physics other than three-particle 
correlation) become the solid data points in Fig. Q] As seen, both the LS and the US three-particle correlators are 
negative and seem to follow a regular trend. Those remaining effects can be due to two-particle correlation from 
clusters together with cluster anisotropies, as well as any other physics. 
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B. Effect of two-particle correlation from clusters 



We now discuss the effect of two-particle correlation. Our approach is to assume that the only remaining effects 
in the correlator measurements are the combined effects of two-particle correlations and cluster anisotropies, and see 
whether the cluster flow parameters extracted from such an approach are reasonable. If they are unreasonable, then 
there may be new physics, such as local strong parity violation. 

After removing three-particle correlation, the remaining three-particle correlator can now be factorized by Eq. (TTJ) 
because particle c is not correlated with particle pair a and /3 through clusters. However, particle c (hydro-like 
particle) can be still correlated with the cluster if clusters are anisotropic. Two-particle correlation from clusters 
together with cluster anisotropy can give non-zero contribution to the measured three-particle correlator. 

We subtract the estimated three-particle correlation effects from the measured US and LS three-particle correlators, 
(cos(a + (3 — 2c))™g as and (cos(a + (3 — 2c))™° as , respectively, as shown in Fig. [TJ What is left in the three-particle 
correlator is two-particle correlation from clusters. Two-particle correlations can exist between a and f3 (and c is not 
from the cluster), or similarly between a and c, or between f3 and c. The former is simply 

(cos(a + (3- 2c)) Q(3 = (cos(A0 Q + Afo + 2(</> clust - rp) - 2(c - i/)))) a p = (cos(A0 a + Ac/>p)) a p ■ v 2 ,ciust • ife.o (40) 

while the latter, e.g. between a and c, is given by 

(cos(a + f3 - 2c)) ac = (cos(A0 Q - 2A0 C - (<£ciust - ip) + (P ~ ip)))ac = (cos(A0 Q - 2A0 C )) QC • ui, c ; us t • vi,p. (41) 

Since direct flow v\ is generally much smaller than elliptic flow vi at mid-rapidity, the latter correlations can be 
neglected and we shall only focus on two-particle correlation effect between a and 0. 

We divide the three-particle correlation corrected results by V2,c used in Ref. [j| to obtain (cos(a + (3 — 2ij))) us and 
(cos(a + j3 — 2tp)) LS . The dilution factors are properly taken into account. Assuming the only remaining correlation 
is from clusters (no new physics), then 

(coa(a + P-2fl>)) ua = x us {cos(a + /3 - 2^))) SA + (1 - x us ){cos(a + - 2i/>)) BB , (42) 

(cos(a + /3-2V)) iS = x LS {cos(a + 0-2i;)) SA + (l-x LS ){cos{a + j3-2i(>)) BB , (43) 

where 

(cos(a + (3 - 2ip)) SA = (cos(A0 Q + A4> p + 20 clust - 2ip)) SA = (cos(a - f3)) SA v 2 SA = w SA v 2 SA , (44) 

(cos(a +j3- 2iJ;)) BB = (cos(A0 Q + A(/)p + 2(f> c i ust - 2<ip)) BB = (cos(a - 0)) BB v 2 BB = -w BB v 2 BB . (45) 

Here v 2 SA and v 2 BB arc the elliptic flow parameters of SA- and BB-clusters weighted by the number of particle 
pairs per cluster (i.e. anisotropy of clusters with each cluster counted by the number of pairs in the cluster). In 
Eqs. (j44|) and (j45|) . we have assumed that the cluster two-particle azimuthal correlation shape is independent of the 
cluster orientation with respect to the reaction plane, so that we can factorize the cluster azimuthal spread and the 
cluster anisotropy. If particle azimuthal distribution in clusters depends on the cluster orientation, then the cluster 
anisotropy should be taken as an effective average anisotropy. From Eqs. (|42|) . p3| . ([44)) and d45l . we have 

(cos(a + /3-2V')) Lrs = x us w SA v 2SA -(l-x us )w BB v 2BB , (46) 
(cos(a + - 2ip)) LS = x LS w SA v 2 SA - (1 - x LS )w BB v 2 BB . (47) 

From Eqs. (j46)) and (|47|) we solve for v 2 SA and v 2 BB : 

1 {l-x L3 ){coB{a + p-2ij,)) us -{l-x us ){cas{a + p-2i ) )) LS 



x LS (cos(a + 13 - 2ip)) us - x us (cos(a + p- 2ip)) 1 



;(1 x ls) x ls(^ x us) 



(48) 
(49) 



Our strategy now is to see what values of v 2 SA and v 2 BB are needed in order to reproduce the correlator measure- 
ments (cos(q!+/3 — 2ip)) us and (cos(a + (3 — 2ij))) LS from STAR and judge whether the required cluster anisotropies 
are reasonable. Table I shows the obtained flow parameters v 2 SA and v 2 BB for SA- and BB-cluster particle pairs, 
respectively. The errors are propagated from statistical errors on the three-particle correlator measurements and 5% 
errors on low-px two-particle angular correlation measurements (same-side Gaussian amplitude, same-side Gaussian 
a, and dipole amplitude). All errors are treated as uncorrelated. 
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TABLE I: Results from cluster model calculation. The x us is the fraction of US pairs from SA-clusters; the x LS is that 
of LS pairs from SA-clusters; the x us and x LS are obtained from (cos(a — j3)) measurements by Eqs. (|22p and (|23p using 
inputs from two-particle angular correlation measurements [8J, assuming no charge difference in BB-clusters (Eq. 1)140 ). Cluster 
(cos(a + P — 2c)) c i ust is the calculated three-particle correlation effect assuming independent particle emission in clusters and 
binary collision scaling for the number of clusters. The v 2 SA and v 2 B „ are the elliptic flow parameters by Eqs. (|48|) and (I49|) 
to reproduce the three-particle correlator results measured by STAR j5|, after cluster three-particle effect removed; errors are 
propagated from statistical errors on the correlator measurements [f| and the assumed 5% error on the two-particle angular 
correlation measurements |g|. The last two columns list the v 2 SA and v 2 BB parameters needed to reproduce the three-particle 
correlator results [§}, assuming vanishing three-particle correlation from clusters. 
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Figure [2] (upper left panel) depicts the obtained v 2SA and v 2BB . The SA-cluster particle pairs i>2 is somewhat 
negative. Negative SA-cluster particle pair v2 is not impossible, perhaps even natural in the jet-quenching picture - 
high-pT particles are suppressed more in the out-of-plane direction, generating more low-px particles, thereby more 
SA-cluster particle pairs out-of-plane than in-plane. Positive anisotropy for BB-cluster particle pairs implies larger 
survival probability of BB-pairs in-plane than out-of-plane, again consistent with the jet-quenching picture. The 
magnitudes of the obtained v 2 SA and v 2 BB seem reasonable, however, the trends towards the most central collisions 
seem unreasonable. More discussions on the extracted flow parameters can be found in Section IIIII 

C. Dependence of results on model assumptions 

We have made two major assumptions in our study as we have italicized in Section [II Al 

(i) The number of clusters scales with binary collisions so we may estimate the cluster size; and 

(ii) Particle emission azimuths within a cluster are independent of each other so we may factorize three-particle 
correlation as the product of two-particle correlations. 

Assumption (i) does not affect the two-particle dilution factor of Eqs. (fl"5j) and (fT6|) because the number of correlated 
particle pairs N c \ ust M^ lust is one of the measured quantities, as we noted already. However, the assumption does have 
an effect on the estimation of the number of correlated triplets N c \ ust M^ lust , thereby on the three-particle dilution 
factor of Eqs. (1551) and A larger A^iiist than the binary scaling estimate would result in a smaller cluster size 

M c \ ust , hence smaller three-particle correlation effect; a smaller A^iiist would result in a larger three-particle effect. 
With binary scaling estimation of iV c i us t, the cluster size varies in the range of -M c i us t ~ 5-10 from peripheral to central 
Au±Au collisions, and the fraction of particles from clusters varies from 5-20% of all particles measured in the final 
state H,|{|- Other studies, from multiplicity correlations, indicate that cluster size is only ~ 3 whereas more particles 
originate from clusters [l3j. If cluster size is Mciust = 3 independent of centrality, then the fraction of particles from 
clusters would be 8-70% of all particles from peripheral to central collisions and the number of clusters would scale 
more strongly than binary collisions. 

Assumption (ii) affects the factorization approximation in deriving the cluster three-particle correlation by Eqs. ([24]) 
and (|23|) . If independent emission of particles from clusters does not hold, then the three-particle cluster correla- 
tion result may not be valid. The factorization approximation also affects the connection between the two-particle 
correlations (cos(A0 Q + A<f>p)) = (cos(A0 Q - A^)) in Eqs. flU) and (H5|) . 

To get a "feeling" about the effects of three-particle correlations from smaller clusters, we repeat our analysis by 
fixing the cluster size to be M c \ us t = 3, independent of centrality, while keeping the measured number of cluster 
particle pairs N c \ ust M^ lust the same. This changes the three-particle correlation effects from clusters. The cluster 
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FIG. 2: (Color online) The v 2 SA and v 2 BB parameters (anisotropies of cluster population weighted by the number of par- 
ticle pairs per cluster) required to reproduce the STAR three-particle correlator measurements [fj after removing cluster 
three-particle correlation effects calculated from two-particle angular correlation measurements || and two-particle correlator 
measurements [jj. The upper left panel assumes binary scaling for the number of clusters; the upper right panel assumes fixed 
cluster size of 3 particles; the lower left panel assumes that cluster size is a factor of 3 smaller than those in the upper left panel; 
and the lower right panel assumes a vanishing three-particle correlation from clusters. The particle elliptic flow measured by 
the event-plane method is shown in the dashed curve, and that from fit to two-particle correlation measurement [H, [12] is 
shown in the solid curve. 



v 2 SA and v 2 BB that are needed to reproduce the measured three-particle correlators change accordingly. Figure [2] 
(upper right panel) shows the obtained v 2 3A and v 2 BB . Both the v 2 SA and v 2 BB increase from the case in the upper 
left panel where binary collision scaling is assumed for the number of clusters resulting in large clusters. Furthermore, 
the values of v 2 SA become mostly positive. The magnitude of v 2 SA appears reasonable. The magnitude of v 2 BB 
seems reasonable for most centralities except peripheral collisions where it is too large. 

Besides the fixed cluster size, we have also tried reducing the cluster size by a constant factor of 3 for all centralities. 
The obtained v 2 SA and v 2 BB (which are needed to explain the three-particle correlator measurements) are shown in 
Fig. [2] (lower left panel). The results are similar to the previous results where cluster size is fixed to 3. 

Reducing the cluster size reduces the effect of cluster three-particle correlation. The extreme would be to assume a 
vanishing three-particle correlation from clusters, and the only contribution to the measured three-particle correlator 
is the combined effect of two-particle correlation from clusters and cluster anisotropy. To test this extreme, we set 
the cluster three-particle correlation to zero, and repeat our analysis to extract the values of v 2 SA and v 2 BB that 
arc needed to fully account for the measured three-particle correlators. The extracted v 2 SA and v 2 BB are tabulated 
in the last two columns of Table Q] and are shown in the lower right panel of Fig. [5] Both the v 2 SA and v 2 BB are 
larger than the other cases, as expected from their trend with reducing effect from cluster three-particle correlation. 
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Both v 2 SA and v 2 BB are positive, and v 2 BB is larger than v 2 SA . The magnitude of v 2 SA seems reasonable (see more 
discussion on this in Section HIT]) . The magnitude of BB seems reasonable in central collisions, but is too large in 
peripheral collisions. 

The signs and relative magnitudes of v 2 SA and v 2 BB in Fig. [5] (lower right panel) can be understood as follows. 
The US three-particle correlator is nearly zero. The three-particle correlators due to SA- and BB-cluster correlations 
should cancel each other. Because particle pairs from SA- and BB-clusters at the same location relative to the reaction 
plane give opposite sign three-particle correlators, the anisotropy of the SA- and BB-clusters have to be of the same 
sign (either positive or negative). Because more US pairs come from SA-clusters than BB-clusters, the anisotropy 
of BB-cluster particle pairs has to be larger than those of SA-clusters in order to have the averaged three-particle 
correlator to be nearly zero. The sign of the LS three-particle correlator is decided by BB-clusters, because more LS 
pairs come from BB-clusters and v 2 BB is larger than v 2 SA in absolute magnitude. BB-cluster particle pairs in-plane 
gives negative (cos(a + /3 — 2ip)) and those out-of-plane gives positive (cos(a + j3 — 2ijj)}. In order for the final averaged 
LS pair correlator (cos(a + (3 — 2ip)) LS to be negative, the v 2BB has to be positive. 



III. DISCUSSION ON CLUSTER PARTICLE PAIR ANISOTROPY 



The cluster anisotropy v 2 SA and v 2 BB obtained above are the magnitudes of modulation of the number of clusters 
weighted by the average number of particle pairs per cluster. One way to gauge whether the obtained i>, SA and BB 
are reasonable or not is to see how large the hydro-like particle i>2.bkgd has to be in order, together with cluster V2, to 
reproduce the observed final state particle w™ 63,8 . Since our clusters are large, having on average 5-10 particles, the 
particles from back-to-back pairs are part of those from SA-clusters. We can therefore consider only SA-cluster pair 
v 2 SA . Note this may not be accurate in peripheral collisions where a BB-cluster of a single back-to-back particle pair 
is not part of a SA-cluster pair, thus the BB v 2 BB has to be also considered. But for our purpose, it is sufficient to 
use only SA v 2 SA to get a "feeling". 

The hydro-like particle «2,bkgd can be obtained from the following two scenarios. 

• Assuming cluster size does not vary with respect to the reaction plane, then cluster v 2 SA translates directly into 
cluster particle v-i by 

particle ^ ( g Q) 



where w 2 SA is the angular spread of cluster particles relative to the cluster axis by Eq. (|28|) . The hydro- like 
particle W2,bkgd is then given by 

(1 - /Kbkgd + fw 2 , SA v 2 SA = « 2 mcas , (51) 
where / is the fraction of particles from clusters. 

Assuming cluster orientation is isotropic (e.g. initial hard scattering products), but cluster size varies from 
in-plane to out-of-plane, then approximately half of the cluster v 2 SA translates into cluster particle V2 by 

v^° = w 2tSA v 2 , SA /2. (52) 
The hydro-like particle i>2,bkgd is then given by 

(1 - /Kbkgd + fw 2 , SA v 2 , SA /2 = v™ as . (53) 

We use the ^2{2d} fit to the two-particle angular correlation data Q as the measured w™ 63,8 , because the non-flow 
same-side correlation peak (part of the cluster pair correlation) should be excluded. In fact, if the fit model [8| used 
to separate elliptic flow and non-flow correlations is accurate, then the fit «2{2d} should be the true elliptic flow 
(correlation related to the reaction plane) [9]. Note that the true elliptic flow may not be necessarily equal to the 
hydro-like elliptic flow, but is the net sum of the hydro-like elliptic flow and the product of cluster correlation and 
cluster elliptic flow [9], as shown in Eqs. (jBTj) and (f53|) . 

Figure [3] shows the obtained hydro-like particle «2,bkgd from Eqs. (|51l) and (|53| together with the «2{2d}. The 
«2{ep} measured by the event-plane method, which contains a large contribution from non-flow, is also shown. Its 
difference from i>2{2d} gives a good estimate of uncertainty of all available elliptic flow measurements. The left panel 
of Fig. [3] shows the case where the number of clusters scales with binary collisions corresponding to the upper left 
panel of Fig. [2] The obtained hydro-like particle V2 is larger than the measured particle V2 to account for the negative 
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FIG. 3: (Color online) hydro-like particle particle U2,bk g d that is needed in order, together with the estimated v 2 SA , to reproduce 
the fit i>2{2d}. The left panel corresponds to n, SA in Fig. [2] upper left panel where the number of clusters is assumed to scale 
with binary collisions. The right panel corresponds to v 2 SA in Fig. [5] lower right panel where the three-particle correlation 
from clusters is set to zero. The fraction of cluster particles is kept as same as that for the left panel. 



cluster v 2 SA , but not much larger. The W2,bk g d values seem reasonable, suggesting that the v 2 SA and v 2 BB may be 
reasonable too. 

The right panel of Fig.[3jshows the case where the three-particle correlation from clusters is set to zero, corresponding 
to the lower right panel of Fig. [21 Although the fraction of particles from clusters does not matter for this case because 
the three-particle correlation from clusters is set to zero, we need the fraction to obtain the hydro-like particle t>2,bkgd 
by Eqs. ([51]) and ([55]) . We use the same fraction of particles from clusters as in the above case (assuming binary 
collision scaling for the number of clusters). The calculated hydro-like particle i>2,bkgd is smaller than the measured 
particle v 2 for cluster size independent of the reaction plane to offset the larger cluster particle v^ a ^ lcle (by Eq. ([501) ). 
For isotropic clusters, the hydro-like particle i>2,bkgd is not much different from the measured particle v 2 because the 
cluster particle wP"| lcle due to cluster anisotropy (by Eq. ([52])) approximately equal to the measured particle The 
hydro-like particle ^2,bkgd is reasonable for this case too, again suggesting that the v 2 SA and v 2 BB required to explain 
the three-particle correlator measurements are reasonable. 



IV. CONCLUSIONS 

Cluster model parameters are calculated by Eqs. J7]), ©, and (f^5|) using two-particle angular correlation data Q 
and two-particle azimuth correlator (cos(ev — /?)) measurements from STAR [5j. Fractions of unlike-sign (US) and 
like-sign (LS) particle pairs from small-angle (SA) and back-to-back (BB) clusters are obtained by Eqs. (f22|) and 
([2"5]) assuming no charge difference in BB two-particle correlation. Three-particle correlation effects are estimated 
by Eqs. (f55]) and ([551) assuming independent emission of particles in clusters, and are positive for SA-clusters and 
negative for BB-clusters. The estimated three-particle effects are removed from the three-particle azimuth correlator 
(cos(cv + /3 — 2c)) measurements Q. The remaining correlator magnitudes of (cos(cv 4-/3 — 2-0)) = (cos(cv4-/3 — 2c))/v2, c 
are assumed to come entirely from cluster two-particle correlations (i.e. no new physics), and are used to determine 
elliptic flow parameters of SA- and BB-clusters by Eqs. (H51) and (05]). 

Cluster size is not measured. A wide range of assumptions are made, ranging from binary collision scaling of cluster 
abundance, resulting in large cluster three-particle correlation, to zero cluster size yielding vanishing three-particle 
correlation effect. These assumptions do not affect the cluster two-particle correlation which is constrained by the 
two-particle angular correlation measurements @. The magnitudes of the obtained cluster anisotropy (azimuthal 
modulation in the number of clusters weighted by the number of particle pairs per cluster) to fully account for the 
three-particle correlator measurements [j| seem reasonable, except for peripheral collisions. The hydro-like particle 
flow magnitude, in order to make up to the measured inclusive particle flow together with the cluster particle flow, 
appears reasonable too. 

Cluster particle correlations may originate from (semi-)hard scatterings, ft is therefore natural to expect that cluster 
effect would increase with pt, although the pT-dependence is not studied in this paper due to the lack of py-dependent 
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measurements of two-particle angular correlations. It is worth to note that the measured azimuth correlators indeed 
increase with px [E| > & n observation that cannot be explained by local strong parity violation but may be expected 
from cluster particle correlations. 

In conclusion, our results from conventional physics of cluster particle correlations suggest that no new physics 
is required to explain the two- and three-particle azimuth correlator measurements by STAR. Our conclusion is 
complementary to the experimental findings in Ref. |5j that the azimuth correlator data in itself does not allow to 
conclude on the local strong parity violation. 
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